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Abstract 

Open charm production in 77 collisions is studied with data collected at e + e~ 
centre-of-mass energies from 189 GeV to 202 GeV corresponding to a total inte- 
grated luminosity of 410 pb _1 . The charm cross section (7(77 — > ccX) is measured 
for the first time as a function of the two-photon centre-of-mass energy in the in- 
terval from 5 GeV to 70 GeV and is compared to NLO QCD calculations. 
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1 Introduction 



The cross section for the interaction of two quasi-real photons is described by three components. 
The first "soft" component is described by the Vector Dominance Model ( VDM) , parametrized 
in hadronic phenomenology by Regge poles. A second "direct" component, the point like 
reaction 77 — > qq, can be calculated in QED. Finally, there is a "hard component" (resolved or 
anomalous QCD component) which requires knowledge of the quark and gluon parton density 
functions of the photon. The direct and hard contributions to two-photon interactions can 
be measured in open heavy flavour production. In this case a hard physical scale is given 
by the c- or b-quark mass. At LEP energies, the direct and single resolved processes, shown 
in Figure ^ are predicted to give comparable contributions to the charm production cross 
section pQ, whereas at low energies the direct process dominates. The main contribution to 
the resolved photon cross section is the photon-gluon fusion process 7g — > cc. Contributions to 
charm production arising from VDM processes and from doubly resolved processes are expected 
to be small pQ. 

This letter presents the measurement of the 77 — ► ccX cross section as a function of the 
two-photon centre-of-mass energy W jry in the interval 5 GeV < < 70 GeV. The data 

correspond to a total integrated luminosity C = 410 pb -1 , collected with the L3 detector [2] 
at centre-of-mass energies y/s = 189 — 202 GeV. The inclusive charm production cross section 
o"(e + e~ — > e + e~ccX) was measured by_L3 at y/s = 91 — 202 GeV 0111 ■ Charm quarks are 
identified by their decay into electrono jlj. 

2 Monte Carlo 

The PYTHIA [5, Monte Carlo is used to model the e + e~ — > e + e~hadrons processes. Quarks 
other than b are taken as massless in the corresponding matrix elements The resolved 
process uses the SaSld photon structure function [Zj. The two-photon luminosity function is 
implemented in the equivalent photon approximation (EPA) jH] with a cutoff Q 2 < m 2 , where 
m p is the mass of the p meson. 

The background sources are e + e~ — > e + e~r + r~, e + e~ — > qq, e + e~ — > t + t~ and e + e~ — > 
W+W". These processes are generated by JAMVG |3j, PYTHIA, KORALZ PUJ and KO- 
RALW ^T] respectively. The detector simulation is performed using the GEANT [12J and 
GHEISHA |TH] packages. The Monte Carlo events are reconstructed in the same way as the 
data. Time dependent detector inefficiencies, as monitored during the data taking period, are 
also simulated. 

3 Measurement of cr(77 — ► ccX) 

The measurement of the (7(77 — > ccX) cross section is performed on the high statistics electron 
sample of 2434 events discussed in Reference The backgrounds from annihilation processes 
and two-photon production of tau pairs are estimated to be 0.75% and are subtracted from 
the data. The background from 77 — > bbX events is modelled with the PYTHIA Monte Carlo 
assuming equal contributions from direct and resolved processes. It is also subtracted from the 
data assuming our measured cross section j3]. The charm purity after background subtraction 
is 75%. 

Electron stands for electron or positron throughout this paper. 



2 



The analysis is restricted to events with visible mass W V - 1S > 3 GeV. W V - 1S is calculated 
from the four-momentum vectors of the measured particles, tracks and calorimetric clusters 
including those from the small angle luminosity monitor. The W v i S distribution is corrected for 
trigger efficiency using data with a set of independent triggers. It varies from 94% at W v - m = 3 
GeV to 98% at W vis > 30 GeV. Table [T] shows the selected data events in bins of W vis and 
the different contributions of the signal and the background predicted by the PYTHIA Monte 
Carlo. A comparison of the visible mass distribution of the data after final selection with the 
expectations of the PYTHIA Monte Carlo at = 189 — 202 GeV is shown in Figure El 



3.1 Unfolding procedure 

The hadronic final state is not always fully contained in the detector acceptance. An unfolding 
procedure is hence applied to obtain the 77 centre-of-mass energy W 77 from W v i B . Figure El 
shows the correlation of the W vis average value, (Wvis), with W 77 as predicted by the Monte 
Carlo. The unfolding corrects for both missing particles and detector resolution by the relation: 

N(W^(i)) = J2A 1 N(W vis (j)), (1) 

j 

where iV(W 77 (i)) is the content of the i-th interval of W jry listed in Table El and N(W Y i B (j)) 
is the content of the j-th interval of W V { S listed in Table The matrix A^j is constructed by 
considering for each Monte Carlo event its reconstructed W vis and generated IV 77 [Hj: 

P(W vis (j)\W^))P(W^)) 

13 EiP(w Vls (j)\w 77 (i))P(w JJ (i)V { ) 

where P{W V \ S \W 71 ) is the likelihood of observing the measured value Wvis given a generated 
value of W 77 and P(W 77 ) is the generated W 77 distribution after the selection cuts. 

The unfolding matrix is determined using charm events generated with PYTHIA. The un- 
folding uncertainty is estimated with PYTHIA using five flavour events instead of charm events 
and using the PHOJET J3| Monte Carlo generator with all quark flavours. This comparison 
leads to an estimated uncertainty of 5%. Charm production in the PHOJET generator is im- 
plemented with a high p t threshold which prevents the use of PHOJET for charm efficiency 
studies. 

After unfolding, the events are corrected for efficiency using the ratio between selected and 
generated charm events in each W 77 interval. The events with 3 GeV < W vis < 5 GeV and 
W vis > 70 GeV are used for unfolding, but are excluded from the measurement due to their 
large correction factors and unfolding uncertainty. 



3.2 Cross section measurement 

The cross section Acr(e + e~ — > e + e~ccX) is measured from the number of events corrected for the 
efficiency with the PYTHIA Monte Carlo in each H / 77 interval and the integrated luminosity. 

To extract the cross section 17(77 > ccX) of two real photons, the photon flux £ 77 [H] is 
calculated and the hadronic two-photon process is extrapolated to Q 2 = 0, where Q 2 is the 
virtuality of a photon. This is done through the following relation: 

r dQ 2 dQ 2 

Acx(e + e" -> e + e~ccX) = / 7^ 75/^77 (W^, Qj, Q 2 2 )a yy ^ C c X (W^)dW^. 
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For each W 71 interval a numerical integration is performed over its width and the unmeasured 
Q 2 of the scattered electrons. 

Table El gives the efficiencies and the cross sections Acx(e + e _ — > e + e~ccX) and (7(77 — > ccX) 
as a function of W 7J . The systematic uncertainties arise mainly from the charm purity estimate, 
followed by the unfolding procedure, the charm efficiency, the beauty cross section [3], the 
photon flux and the trigger efficiency, all summarized in Table 01 The photon flux uncertainty 
is estimated by the comparison of the two models (H1IIE]. The uncertainty due to the other 
background processes is negligible. Due to the unfolding, the data points are correlated. The 
correlation matrix is given in Table 



4 Comparison with theory and interpretation 

Figure |U shows the differential cross section A<r(e + e - — > e + e~ccX)/Al / l / 77 as a function of W 77 
measured in the interval 5 GeV < TU 77 < 70 GeV. The expected slope from the PYTHIA Monte 
Carlo is steeper than that of the data. 

Figure El compares the measured cross section (1(77 — > ccX) as a function of W 77 with NLO 
QCD calculations . The calculations use massive quarks in the matrix elements. The charm 
mass, m c , is fixed to 1.2 GeV, the renormalization and factorization scales are set to m c and 
2m c , respectively, the QCD parameter A^ CD is set at 227.5 MeV, and the GRS-HO [TH] photon 
parton density function is used. Using this set of input parameters, the NLO QCD predictions 
reproduce well the energy dependence and the normalization. The calculation with m c = 1.5 
GeV results in about 50% lower cross section values, except the first point, where it is lower by 
25%. A change in the renormalization scale from m c to 2m c decreases the QCD prediction by 
10% and 30% at low and high W 77 respectively. 

We compare also the measured charm cross section with the total cross section of hadron 
production in two-photon collisions [T^], scaled by an arbitrary factor 1/20. The slope of 
(7(77 — y ccX) is clearly larger than that of 17(77 — > hadrons). 

A more quantitative comparison results from fits to the data. A parametrisation [20] of 
the form o = As e + Bs~ v , "Pomeron + Reggeon", with s = 1V 77 describes well the energy 
behaviour of all the total hadron-hadron cross sections with universal values e = 0.093 ± 0.002 
and 1] = 0.358 ± 0.015 JH]. A fit to our data with e, A and B being free parameters and with 
fixed r] = 0.358 yields in the interval 5 GeV < W 77 < 70 GeV: 

e = 0.40 ± 0.03 (stat.) ± 0.07 (syst.) 

A = 1.3 ±0.3 (stat.) ± 0.7 (syst.) nb 

B =44.0 ±3.8 (stat.) ±11.1 (syst.) nb 

X 2 /d.o.f. = 8.9/2 

Correlations between the data points are taken into account. The systematic uncertainty on 
each parameter is estimated by performing fits with total uncertainties and statistical uncertain- 
ties only in \ 2 - The systematic uncertainties are then evaluated by subtracting in quadrature 
the uncertainty of the statistical only fit from that given by the total uncertainties. The fitted 
value of e is higher than the universal value [To] . 

Modifications to the parametrization of the total cross section were recently proposed |21j . 
In order to allow for a larger slope for the high energy behaviour, an admixture of "hard" and 
"soft" Pomerons with the Pomeron powers e of 0.4 and 0.1 respectively was introduced. Our 
data correspond to a direct measurement of the "hard" component of photon-photon collisions. 
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A fit to the data by the form a = A s e , performed in the interval 10 GeV < < 70 GeV, 
yields: 



The value of e is in agreement with a similar fit to the measurements of the photoproduction 
of J mesons at HERA [221123]. 
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W vis 


(GeV) 




iVcc 


AU 




-Wbkgd 


3 


- 5 


261 


208 


48 


0.2 


4.6 


5 


- 10 


741 


519 


186 


31 


5.3 


10 


- 15 


527 


360 


116 


50 


1.3 


15 


-25 


541 


352 


127 


61 


1.1 


25 


-40 


276 


183 


55 


36 


2.3 


40 


- 70 


85 


47 


22 


12 


3.8 




> 70 


3 


0.3 


2 


0.5 


0.2 



Table 1: The number of events for the data, N data , the expected signal, iV cg , the background 
from light flavours, iV u ds, b production iVb and other processes iVbkgd as a function of W v i S . 



A1U 77 ( 


GeV) 


N C c unfolded 


Efficiency (%) 


Aa(e+e- -> e+e^ccX) (nb) 


a (77 -> ccX) (nb) 


5 - 


10 


428 


0.33 ± 0.01 


0.316 ± 0.013 ± 0.040 


21.0 ± 0.8 ± 2.9 


10- 


15 


268 


0.56 ± 0.02 


0.117 ± 0.004 ± 0.015 


18.6 ± 0.7 ± 2.6 


15- 


25 


389 


0.70 ± 0.02 


0.136 ± 0.005 ± 0.018 


24.2 ± 0.9 ± 3.4 


25- 


40 


325 


0.72 ± 0.03 


0.110 ± 0.004 ± 0.015 


33.0 ± 1.3 ± 4.9 


40- 


70 


196 


0.57 ± 0.03 


0.084 ± 0.004 ± 0.013 


38.8 ± 1.9 ± 6.2 



Table 2: Unfolded number of charm events, iV cg , efficiencies and cross section values as a 
function of W 77 . The first uncertainty on the cross section is statistical and the second is 
systematic. The cross section values are evaluated at the centre of each W 77 interval. 





1U 77 


Source of uncertainty 


5-10 GeV 


10 - 15 GeV 


15 - 25 GeV 


25 - 40 GeV 


40 - 70 GeV 


Charm purity 


11.3 


12.2 


12.2 


12.6 


12.6 


Unfolding 


5.0 


5.0 


5.0 


5.0 


5.0 


Charm efficiency 


3.8 


3.7 


3.3 


3.9 


5.0 


a(e+e- -> e+e-bbX) 


3.5 


1.6 


2.8 


3.4 


4.1 


Photon flux 


0.1 


1.0 


1.5 


2.7 


4.8 


Trigger efficiency 


2.1 


2.1 


2.2 


2.3 


2.3 


Total 


13.6 


14.0 


14.1 


14.9 


15.9 



Table 3: Systematic uncertainties (in percent) on the cross section for the process 77 — > ccX. 



W, n (GeV) 


5-10 


10 - 15 


15-25 


25-40 


40-70 


5-10 


1.0 










10- 15 


0.739 


1.0 








15-25 


0.294 


0.703 


1.0 






25-40 


0.090 


0.302 


0.757 


1.0 




40-70 


0.014 


0.093 


0.388 


0.778 


1.0 



Table 4: Correlation matrix of the data after unfolding. 
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spectator 
jet 



Direct Single Resolved 

Figure 1: Diagrams contributing to charm production in 77 collisions. 
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W vls (GeV) 

Figure 2: The visible mass spectrum, W v i s , for the inclusive electron data at y/s = 189—202 GeV 
compared to the PYTHIA prediction. The Monte Carlo spectrum with all flavour contributions 
is normalized to the integrated luminosity of the data after scaling to the measured charm and 
beauty cross sections according to our measurements [3]. Other background sources are also 
shown. 
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W w (GeV) 

Figure 3: Profile of the distribution of W^is as a function of W 77 using the PYTHIA Monte 
Carlo. The open circles correspond to the mean values and the error bars represent the r.m.s. 
of the distribution. The dashed line corresponds to perfect correlation. 
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W w (GeV) 



Figure 4: The cross section Aa(e + e — > e + e ccX) / AW 77 as a function of at y/s = 189—202 
GeV. The dashed line corresponds to the leading order PYTHIA Monte Carlo prediction. 
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Figure 5: The cross section (1(77 — > ccX) as a function of W 77 at ^/s = 189 — 202 GeV. 
The dotted curve is the total cross section 17(77 — > hadrons) measured by L3 [TH] scaled by an 
arbitrary factor 1/20. The continuous line is the NLO QCD prediction, while the dashed-dotted 
and dashed curves show the expectation from the direct and resolved process respectively. 
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